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IR finite one— loop box scalar integral with massless internal lines 
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Theoretical Physics Division, Rudjer Boskovic Institute, P.O. Box 180, HR-10002 Zagreb, Croatia 

The IR finite one-loop box scalar integral with massless internal lines has been recalculated. The 
result is very compact, simple and valid for arbitrary values of the relevant kinematic variables. It 
is given in terms of only two dilogarithms and a few logarithms, all of very simple arguments. 
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I. INTRODUCTION 

Exclusive hadronic processes at large momentum 
transfer in which the total of particles (partons) in the 
initial and final states is iV > 5 are becoming increas- 
ingly important in testing various aspects of QCD. The 
tree-level results for the cross section of these processes, 
being proportional to a high power of the QCD coupling 
constant ctg, are very sensitive to the variation of the 
renormalization scale and scheme being used. Conse- 
quently, to stabilize the leading-order (tree-level) predic- 
tions and to achieve a complete confrontation between 
theoretical predictions and experimental data, at least 
one-loop (NLO) corrections are necessary. 

The main theoretical difficulty in obtaining the NLO 
predictions consists in the treatment of the occurring TV- 
point one-loop scalar and tensor integrals with massless 
internal lines. Through the scalarization and reduction 
procedures [I], the computation of such integrals 

reduces to the calculation of a set comprised of six box 
(N = 4) diagrams, one of which is finite while the other 
five are IR divergent. Being of fundamental importance 
for one-loop calculations in pQCD with massless quarks, 
it is very important for practical purposes that the re- 
sults for these integrals are obtained in forms as compact 
as possible (expressed in terms of as fewer as possible 
dilogarithmic functions) with as simple as possible argu- 
ments. 

The IR divergent basic box integrals have been con- 
sidered in Refs. j|, |[ |). A more general and detailed 
analysis has been given in Ref. M. As for the IR finite 
box integral, it has been calculated in Ref. J8fl, and later 
obtained in a more compact form in Ref. jfj . From the 
practical point of view, in addition to not being most 
compact, a disadvantage of this result is a very compli- 
cated structure of the arguments of the occurring dilog- 
arithmic and logarithmic functions. An alternative ex- 
pression for the IR finite box integral has been obtained 
in Ref. [Q, where it has been related to the IR finite 
triangle scalar integral. The result of Ref. expressed 
in terms of only two dilogarithms, however, is not valid 
for arbitrary values of the relevant kinematic variables. 

Making use of the method of Ref. (To), in this work 
we carefully recalculate the IR finite box scalar integral. 
Our result is as simple and compact as that of Ref. ficfl 
and, as the result of Refs. [|[ |) is quite general, i.e. ap- 
plicable to all kinematic regions. In addition, to be more 
compact and general, the advantage of our result over 
the results previously obtained is also the fact that one 



can very easily separate the real and imaginary parts of 
logarithms and dilogarithms, making it more appropriate 
for numerical calculations. 

The paper is organized as follows. Sect. 2 is devoted 
to introducing the notation and to some preliminary con- 
siderations. In Sect. 3, using the Feynman parameter 
method and the Mellin— Barnes integral representations, 
we evaluate the IR finite box scalar integral by relating it 
to the IR finite triangle scalar integral. In Sect. 4 we give 
some concluding remarks. The details of the evaluation 
of the IR finite triangle scalar integral are given in the 
Appendix. 



II. PRELIMINARIES 

The scalar one— loop box integral with massless inter- 
nal lines in D— dimensional space— time is given by 



h(pi,P2,P3,Pi) = (M ) 



d D l 



1 



{2tt) d A 1 A 2 A 3 A 4 



(1) 

where pi (i=l,2,3,4) are the external momenta, I is the 
loop momentum and ii is the usual dimensional regular- 
ization scale. As indicated in Fig. |], all external mo- 
menta are taken to be incoming, so that the massless 
propagators are 



A 2 
A 3 
Ai 



I 2 + ie , 
(l+ Pl ) 2 +ie, 
(I+P1+P2) 2 +ie, 
(I +P1+P2 +P?,) 2 + ie • 



(2) 



where the quantity ie (e > 0), represents an infinitesimal 
imaginary part specifying on which side of the cut a mul- 
tivalued function has to be evaluated. We take the cut 
along the negative real axis. 

Combining the denominators using the Feynman pa- 
rametrization formula, performing the D— dimensional 
loop momentum integration, introducing the external 
"masses" p 2 — mf (i = 1,2,3,4) and the Mandclstam 
variables s = (pi + P2) 2 and t = (p2 +P3) 2 , the integral 
(|l|) becomes 



t r 4. 2 2 2 2\ 



i r(4-D/2) 



(4tt) 2 (4n^i 2 ) D / 2 - 2 

/ dxida^dxsda^ 8{x\ + X2 + x 3 + £4 — 1) 
Jo 



2 



Pi 



P2 



l+Pl 




1+P1+P2 



Pa 



I+P1+P2+ P3 



FIG. 1: One-loop box diagram. 
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Depending on the number of external massless lines, 
there are six special cases of the above integral. Adopting 
the notation of Ref. , we denote them by 



rflm 
1 4 


= 7 4 ( S ,i;0,0,0,0,), 


(4) 


jlm 
1 4 


= h(s,t; 0,0,0, ml), 


(5) 


j2me 
1 4 


= h(s,t; 0,7712,0,7714), 


(6) 


j2mh 
1 4 


= h(s,t; 0,0, 7773,7774), 


(7) 


r3m 
M 


= 7 4 (s,i; 0,7772, 7773, 7774), 


(8) 


j4m 
M 


r / , 2 2 2 2\ 
= -/ 4 (s, 7; m 1; 777 2 , 777 3 , 777 4 ). 


(9) 



The above integrals constitute the fundamental set of in- 
tegrals, in the sense that an arbitrary one— loop N(> 
5)— point integral with massless internal lines can be 
uniquely expressed as a linear combination of these in- 
tegrals with the coefficients being rational functions of 
the relevant kinematic variables. The integrals 1|) are 
IR divergent and need to be evaluated in D = 4 + 2em 
(sir > 0) dimensions 0. On the other hand, the in- 
tegral (||) is IR finite and can be calculated in D = 4 
dimensions. 

Setting D = 4 in (0), we obtain 



7/4- 2 2 2 „,2\ 

l4\S, t, 777 l 7 777 2 , 777 3 , 777 4 ) 



da;i dx2 da;3 da:4 



(47T)" JO 

x 5{x\ + X2 + X3 + X4 — 1) (xia;3 s + X2X4 t 

+ X\X2 m\ + X2X3 m\ + X3X4 7773 + 2:1X4 7774 + ie) 2 . 

(10) 

This integral is invariant under some permutations of ex- 
ternal "masses" and Mandelstam variables. Namely, in- 
troducing the following set of ordered pairs 



{ (s,t), (m\,m%), (7772,7774) }, 



(11) 



one can easily see that the integral (|10|) is invariant under 
the permutations of ordered pairs given in (|ll"l), as well as 
under the simultaneous exchange of places of elements in 
any two pairs. Furthermore, the integral is invariant un- 
der the simultaneous change of the signs of the kinematic 
variables and the causal ie. Consequently, it is sufficient 
to analyze only the cases of the integral ( |l0| ) in which the 
number of positive kinematic variables is larger than the 
number of negative ones. 



III. CALCULATION AND RESULTS 

Eliminating the S— function in ( |l(i| ) by performing the 
x\ integration, we get 



T I j. 2 2 2 2\ 
7 4 (S, t; 777i, 777 2 , 777 3 , 777 4 ) 



(Wo Vo 



1 — X\ pl~Xi—X2 



dxi / d2;2 



dxv 



x [(1 - 2:1 - 2; 2 - 2:3) (2:1 7774 + 2; 2 t + x 3 7773 ) 

+ 2:12:2 777i + 2:i.T3 S + 2:22:3 7772 + ie] ' • (12) 

Next, if we first make the substitution 2:3 — * (1 — x± — 
2:2)2:3 and then 2:2 — > (1 — %i)%2> we nn( l that 



T ( + 2 2 2 2^ 

i4(s, t; m\, m 2 , m 3 , m 4 ) 



(4tt) 2 



1 ,1 ,1 

dx\ / dx2 / d2;3 
[i Jo Jo 



x (1 - £2) j(l - xi)(l - x 2 ) [x 2 x 3 ml + 2:2(1 - 2:3) t 
+ (1 - 2:2)2:3(1 - 2:3) 7773] + xi [x 2 mf + (1 - 2:2)2:3 s 



+ {l-x 2 )(l-x 3 )ml] +ie 



(13) 



Performing the x\ integration and a few simple rear- 
rangements, the integral takes the form 



WJo dX2 Jo dX3 



T I + 2 2 2 2\ 
i 4 (S, i; 777i, 777 2 , 777 3 , 777 4 ) 

x [x 2 m\ + (1 - 2:2)2:3 s + (1 - 2; 2 )(1 - 2:3) 777| + ie] 

x [2:22:3 "72 + £2(1 - £3) t+(l- 2:2)2:3(1 - 2:3) 7773 + ie] 1 

(14) 

Next, doing the 2:2 integration and expressing the remain- 
ing integral in terms of logarithms and dilogarithms, one 
arrives at the result obtained in Ref. |)|, which is valid 
for all values of the kinematic variables. The result of 
Ref. Q is given in terms of logarithmic and dilogarith- 
mic functions with the argument of the form 



1 + (h ± V/4 + i/se 

J2 - ie V 



where ft (i ~ 1,2,3,4,5) are real functions of the kine- 
matic variables. As seen from the above, the structure of 
the argument is very complicated since ie appears both 
outside and inside the square root. There are two con- 
sequences of this: first, it is very dificult to determine 
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the sign of the infinitesimal imaginary part of the ar- 
gument; second, it is extremely complicated to use the 
dilogarithmic identities with the aim of simplifying the 
final result. It should also be pointed out that since the 
final result contains products of logarithms, the signs of 
the imaginary parts of arguments are important not only 
for determining the imaginary part of the result, but also 
for the real part of the result. From the practical point 
of view, this represents a disadvantage of the result for 
the IR finite massless box scalar integral derived in Ref. 

To avoid the above-mentioned problems, we take a dif- 
ferent approach and proceed to evaluate the box integral 
by expressing it using the Mellin-Barnes integral repre- 
sentation To this end, we first need to find the 
Mellin-Barnes representation of the following expression: 

J = ( VW2 a + yi(l -y 2 )b + (1 - yi)y 2 (l -2/2) c + ie) -1 , 

where y\, y 2 G [0,1], and a, b and c are arbitrary real 
variables different from zero. In what follows we also 
need the Mellin-Barnes representation of the expression 

{l-z)~ k 



/>7+ioc z^'+ioo 

x / ds ds' r(s)iV) r(i - s - s') 

J -y — ioo J V — ioo 



(l-z)- 



1 



7+100 



dsT(s)T(k - s){-z)- 



2niT(k) J7 _ ioa 
0<7<Refc, |arg(-z)| <tt. (16) 



Rewriting the expression ( |15| ) as 
J 



yi[y2a+ (1 - y 2 )b + ie] 

xi "i+ j-jg wa-^c+i^ s (17) 

yi yi a + (1 - 2/2) + ie 



and utilizing the relation (|16|), J can be written in the 
form 



2 /■7+ico 

J = - — - / ds 



2m 



r(«)r(i-*) y ;- 



7—100 



[(i-yi)(y 2 (i-y 2 )c + ie)] s 



x[y 2 a + (1 -y 2 )b + ie] 



< 7 < 1. (18) 



Now, applying the relation ([T^) to the factor (y 2 a + (1 — 
2/2) & + ie) s ~\ appearing on the right-hand side of (18), 
one finds that 

1 1 i-i+ico i-f'+ico 

J = 77T^ — - / ds / ds'r( s )r(s')r(i-s-s') 

(2tti) 2 a + ie L_ ioo L_ ioo 



c + ie 



b + ie 



(l- yi ) s (l-y 2 ) s + s ' \a + iej V a + ie , 

0<7, 7 ' <1, 7 + 7' < 1, (19) 

which represents the desired Mellin-Barnes representa- 
tion for J given in(15f). If use is made of the relation 
(|l9|), the integral (|l4 ) becomes 



j / 1 2 2 2 2^ 
i4(s, t; m 1 , m 2 , to 3 , to 4 ) 



1 



1 



(4tt) 2 (27ri) 2 m| 



7 — 100 100 
2 



mi+ieV" / t + ie ^ s rl 



in; + le / \ m; -1- le 



Jo 



dx 2 / dx 3 



•'-2 ^ <■,> ' ' 

x 2 to 2 + (1 - x 2 )x 3 s + (1 - x 2 )(l - x 3 ) to| + ie' 

(20) 

After performing the x 2 and x 3 integrations with the help 
of the formula 



r(a)(l-z) b 
Rea>Refe>0, |arg(l - z)| < tt, (21) 



one arrives at 



t / 1 2 2 „,2 „2\ 

i4(s, t; m l7 to 2 , to 3 , to 4 ) 



i 1 



1 



(4tt) 2 (27ri) 2 (ml + ie)(m 2 +ie) 



/•7-t-ioo />7 -1-100 

x / ds ds' [r(s)r(s')r(i - s - s ')f 

J 7 — ioo J 7'— ioo 

(to 2 + ie)- s (to 2 + ie)~ s (s + ie) _s ' (t + ie)" 



(m| + ie)~ s (m| + ie)~ s (ml + ie)~ s ' (m| + ie)~ s ' ' 
0< 7 , 7 '<1, 7 + 7' <1, (22) 

which is the sought Mellin-Barnes representation of the 
IR finite box scalar integral. 

Let us now make a digression and, for a moment, con- 
sider the IR finite scalar one-loop triangle integral with 
massless internal lines in 4-dimensional space-time. It 
corresponds to the Feynman diagram of Fig. |^, and is 
given by 



^3(91,92,93) 
d 4 k 



(23) 



(2iry [k 2 + ie] [(k + 91) 2 + ie] [(k + 91 + q 2 ? + ie] ' 

where g,-, i=l,2,3 are the incoming external momenta. 
Upon combining the denominators with the help of the 
Feynman parametrization formula, integrating out the 
loop momentum and introducing the external "masses" 
9 2 = fj, 2 (i = 1, 2, 3), the integral (|2^ ) takes the form 

h(t4, m!> M§) = 1 w / dxi dx 2 dx 3 S(xi + x 2 + x 3 - 1) 

x (X1X2 /ii + X2X3 /x 2 + X1X3 /i 3 + ie) -1 . (24) 

In this form it is evident that the triangle integral is in- 
variant under arbitrary permutations of external masses 

After eliminating the <5-function by performing the x 3 
integration, introducing the new variable x 2 — > (1— X\)x% 
and taking into account Eq. (|l^), one can write 



1 



1 



(4tt) 2 (27ri) 2 /i 2 + ie 



7+100 p-y +100 

ds / ds' 

7— ioo J 7' — ioo 
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FIG. 2: One-loop triangle diagram. 



1C 



xT(s)T(s')T(l - s - s' 

x / dx^l^il - x^ 3 / dx 2 x S 2~ 1 (l - x 2 ) 
'a Jo 



(25) 



Next, performing the x\ and x 2 integrations with the 
help of the formula (21) (with z — 0), one obtains 



iM,i£,i&) 



i 



(4tt) 2 (27ri) 2 /x 2 + ie 



7+ioo /«7 +ioo 

ds / ds' 



7—100 



7—100 



x[r( s )r( s ')r(i- s - s ')] 2( ^ + ier(/i ' + ie) 



(rf+ie)- s (,i 2 +ie)-*" 
0<7,y<l, 7 + 7' < 1, (26) 

which represents the Mellin-Barnes representation of the 
triangle scalar integral. 

If we now compare the Mellin-Barnes representations 
of the scalar box and triangle integrals given by Eqs. J22]) 
and (pfj), respectively, we arrive at the relation 



I±(s, t; mf, m 2 , ra 3 , m 4 ) = 7 3 (si, m 2 m 3 , to 2 .™!): (27) 

which is valid provided the following condition is satis- 
fied: 



(it + ie) s (v + ie)~ 



(uv + ie) 



(28) 



where (u, v) is any element of the set of ordered pairs 
in ( |TT| ) and s" stands for s or s'. As is easily seen, the 
above condition does not hold only if u and v are both 
negative (u, v < 0). Since, as mentioned earlier, we are 
considering the cases with more positive than negative 
parameters, it is clear that the case when ( |2§| ) is not 
valid can appear at most in one of the pairs. 

Let us now examine in detail the product appearing 
on the left-hand side in (^8|), corresponding to the case 
when u, v < 0. One can then write 



(u + iey s (v + ieY 



exp[— 2m sign(e) s"] 



— (uv) s [cos(7r s") — isign(e) sin(7r s")] 2 
= (uv) ~ s [1 — 2i sign(e) sin(7T s") cos(n s") 

-2sin 2 (7rs")] 
= (uv)~ s {1 — 2i sign(e) sin(7r s") exp[— i7rsign(e) s"]} 
= (uv + iey s " - 2isign(e) sin(vrs") (-uv + ie) _s "(29) 
Next, with the help of the formula 

T(z)r(l -z) = 7rsin _1 (7rz), 

one finds that 

(u + ie)~ s "(u + ie)~ s " = (uv + ie) _s " 
2i7r sign(e) 



— (-uv + ie)- 



(30) 



r(s")r(i - s") 

Owing to the symmetry of the pairs, with no loss of 
generality, we can now suppose that the "problematic 
pair" (i.e. the pair in which both masses are negative) is 
(mf,m|). In this case, the first term on the right-hand 
side in ( pOj ) is of the form that makes it possible to rep- 
resent the box integral in terms of the triangle integral, 
while the second term is a correction originating from the 
fact that the first term does not keep the information re- 
garding the signs of u and v. 

Let us now evaluate the correction term. Upon com- 
bining Eqs. (|22| ) and (|30|), the Mellin-Barnes representa- 
tion of the correction term is 



TSt .222 „,2\ 

K(s,t,m 1 ,m 3 ,m 2 ,m 4 ) = 



i 1 — 2i7rsign(e) 
(4-n-) 2 (27ri) 2 ra\ra\ + ie 



7+ioo 

ds 



7+i ~,r\£) [F( s ')r(i- S - S ')r 



7 — 100 ^7—100 
2,™2 1 • \-s 



T(l-s) 



(-m{ml + ie)~ s (st + ie)~ s 
(m 2 m1 + ie) _s (m^m 2 + ie)~ 



(31) 



which, after making use of the formula (21) (with z = 0, 
a = 1 — s and b = s'), is found to take the form 

1 



K(s, t, mj, m 3 , m 2 , TO4) = 



-2i7r sign(e) 



(47r) 2 (27ri) 2 m\ml + ie 



1 /"7+ioo y>7 +ioo 

x I dy ds ds'T(s)T(s')T(l - s 

?0 J 7 — ioo J 7'— ioo 
,s'-l ( m 2 m 2 



X- 



V 



(— mf m 3 + ie) s (st + ie) 



(1-1/) 



s+s' 



(m 2 mj 



ie 



f 2 1 



nw-(32) 



ie 



If we now compare this expression with the expressions 
@ and @, with yi = 1/2, we get 



TSt j. 2 2 2 2\ 

K (s, t, m\, m 3 , m 2 , m 4 ) 



-2in sign(e) 



(4tt) 2 



(33) 



x / <^y[ym\m\ + (1 - y) st - y(l - y)m\m\ + ie 



Next, after employing the identity 

1 



V m\m\ + (1 - y) st - y(l - y) 



1 



m 2 m 2 + ie 
1 



y-(yi- ie) y - (y 2 + ie) 



(34) 
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where 



and 



Vl 



2/2 



mfm^ — m^ml + st + A 

m 2 m 2 — m\m\ + st — A 
2m 2 m 2 



are the roots of the equation 

y m\m\ + (1 - y) st - y(l - y) m\m\ = 0, 

while 

A = [(st) 2 + (m 2 m 2 ) 2 + (m 2 2 m 2 A ) 2 - 2 stm 2 ml 

— 2 stm^m 2 — 2 m\m^m\rn 21 \ j (35) 

and performing the remaining integration, we arrive at 
the final result for the correction term: 



K(s, t, m 2 , m 2 ,m 2 , m\ 



-2iir sign(e) 



(4tt) 2 A 



X [in (m 2 m\ + st - (m\m\ — A)(l + ie)) 

- in [m\m\ + st - (m\m% + A)(l - ie))] . (36) 

As it is seen, the correction term K is written in the form 
in which the symmetry with respect to the pairs (s, t) and 
(m^m 2 ) is evident. 

Based on the above consideration, the expression for 
the IR finite one-loop scalar box integral can be written 
as 

t / i 2 2 2 2 \ t I + 2 2 2 2 \ 

l4(s,t;m 1 ,m 2 ,m 3 ,m 4 , e) = I^(st 1 m l m :i ,rri 2 m Al e) 

+ K(s,t,m±,m 2 ,m 2 ,m 2 ,e). (37) 

The first term on the right-hand side in ( |37j ) represents 
the expression for the IR finite triangle scalar integral. It 
is evaluated in detail in the Appendix, and is given by 



i 



(4tt) 2 vz(xi 



X2) 



2Li 2 — 

X2 



2Li 2 ( — ) + In [xix 2 + iesign(i/ 3 ) ] 



-X\ J \ —X2 J J J 

where x\ and X2 are the roots of the equation 
x V\ + (1 — x) v-i — ie(1 — x) 1/3 = 0, 

and 

{^1,^2,^3} =~P{f4,l4,l4}> 

where V denotes the permutation of the masses chosen 
in a such a way that the roots x\ and X2 are outside 
the interval [0, 1]. If the signs of all the masses /z 2 are 



the same, 1^3 should be chosen so as to correspond to 
the mass with the smallest absolute value. On the other 
hand, if the signs of the masses are not all the same, 1^3 
ought to be chosen so as to coincide with the mass whose 
sign is opposite to the sign of the other two masses. 

The result (^) is similar to the result given in Ref. 
jnj, but is more convenient for numerical calculations 
since it is easier to separate the real and imaginary part. 
Note that as a consequence of choosing in an appro- 
priate way, the causal ie appears only in one of the loga- 
rithms in d38|). 

The second term on the right-hand side in (^) is the 
correction term. As stated earlier, it accounts for the fact 
that in the product of two kinematic variables we lose 
information about the phases which can lead to missing 
some poles. The expression for the correction term given 
in ( |3q ) corresponds to a special case characterized by 
the fact that the mf,m 2 < 0. This term, however, can 
be written in a general and for practical purposes more 
convenient form: 

K(ai, Pi,a 2 , ^2,0:3, (3 3 ,e) = 
i -2i7rsign(e) ^ 



(4tt) ; 



x In 



(E 



3& 



- [a l f3 l - A)(l + ie) 
(oik + A)(l - ie) 



(39) 



with A given in (|35|) . It should be observed that in prac- 
tice at most one of the terms in the sum in (|39| ) is different 
from zero. 

The above result for the IR finite scalar box integral, 
as stated earlier, has been derived under the assumption 
that the number of positive kinematic variables is larger 
than the number of negative ones. If the opposite is true, 
then owing to the symmetry 



l4(s, t; m 2 , m. 
= h{-s,-t 



2,mg,m 



'3 1 



(40) 



which is obvious from (|1C|) , the corresponding expression 
can simply be obtained from the one given above simply 
by changing the signs of all six kinematic variables as 
well as of the causal e. 

Finally, we have numerically checked our result for the 
IR finite scalar box integral with massless internal lines 
with the corresponding result of Ref. || and found agree- 
ment in all the kinematic regions. 



IV. CONCLUSION 

Using the Mellin-Barnes integral representations, the 
calculation of the IR finite one-loop box scalar integral 
with massless internal lines has been reduced to the cal- 
culation of the IR finite triangle scalar integral with mass- 
less internal lines. The result is presented in a very simple 
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and compact form (only two dilogarithms) and, owing to 
the fact that we have kept the causal ie systematically 
throughout the calculation, is quite general, i.e. valid for 
arbitrary values of the relevant kinematic variables. From 
the practical point of view, the main advantage of our re- 
sult over the results previously obtained is that the argu- 
ments of the occurring logarithms and dilogarithms are 
very simple, which enables one to separate the real and 
imaginary part of the integral very easily. This makes 
our result for the IR finite one-loop scalar box integral 
more appropriate for numerical calculations. 
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Appendix 

In the Appendix we present the details of the calcula- 
tion of the IR massless scalar one-loop triangle integral 
given in (f38|). 

Eliminating the S— function in (p4{) by performing the 
X3 integration, making the change of variables X2 — ► (1 — 
x\)x2 and performing the x\ integration, the integral (124) 
becomes 



(47T) 2 JO 



d:r 



(41) 



In [x fif + (1 - x) ^§ + ie] — In [x(l — x) fj% + ie] 
x n\ + (1 — x) — x(l — x) ^1 

where we have replaced the variable xi by x. The poles 
of the integrand are 



Xl 2 




t4 



f4 



2 2 

i Mi Ms 

2 2 
A*2 



_ (42) 

As can easily be seen from the denominator in (pd|), the 
integration path will not pass through the pole if one of 
the following two cases is realized: first, if the sign of 
/i| is opposite to the signs of [i\ and fi^; second, if all 
the masses are of the same sign and fi 2 has the smallest 
absolute value. Now, owing to the fact that the integral 
-^(Mi! A*ii A*!) is invariant under the permutations of the 
parameters /it 2 , one can easily show that it is always pos- 
sible to permute the parameters /i 2 in a way that the 
poles x\ and X2 are outside the interval [0, 1]. In trying 
to accomplish that, one of the following four possibilities 
can arise: the poles are complex conjugate, both poles 
are in the interval (— oo, 0), both poles are in the interval 
(1, oo) and one pole is in the interval (— oo, 0) and another 
in the interval (1, oo). In the following considerations we 



assume that the parameters /if are already chosen in such 
a way that the above requirement is satisfied. 

Due to the fact that the poles belong to one of the 
above-mentioned cases, it is convenient to rewrite the 
integral (fill ) in terms of x\ and X2- The integral (|4] 
then takes the form 



i 1 f 1 



1 



x In 



1 + x 



1 — Xl — X2 



X\X2 

In [xiX2 + iesign(/i 



(a; - xi)(x - x 2 ) 
In [x(l - x)} 

(43) 



Observe that the arguments of the first two logarithms 
are positive and that there are no poles on the path of the 
integration. Now, by the partial fraction decomposition 
the integral (M3) can be written as 



1 



(4tt) 2 i4(xi-x 9 ) 
x [Gi(x u x 2 ) + G 2 (x 1 ,X2) + G 3 (x 1 ,X2)}, (44) 



where 
Gi(xt,x 2 ) 

G 2 (x 1 ,x 2 ) 



dx — 

x — X\ 



In 



1 + x 



- / dx- 



ln 



X-X 2 



dx- 



1 + x 



1 - Xi - X2 

X1X2 
1 - X\ — X2 



XiX 2 



(,45) 



x — X\ 



•In [x(l -x)\ 



dx- 



1 



/o x - x 2 
G 3 (xi,x 2 ) = In [xix 2 + iesign(^)] 



■ In [x(l - x)] , (46) 



dx- 



1 



X X\ Jq X — X2 



dx- 



(47) 



Let us now evaluate the above integrals in turn. Chang- 
ing the integration variable 



1 + x 



1 — Xl — X2 1 — X\ 



X\X2 



X 2 



in the first, and 



1 + x 



I — X1—X2 l — X 2 



X\X2 Xl 

in the second integral in ([45]), one finds 

j-*3 



Gi(xi,x 2 ) 



dx- 



1 



x- 1 



In 



dx- 



x - 1 



In 



1 — Xl 

x 2 
1 - x 2 



xi 



(48) 
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Next, after passing to the new integration variable x — > 
1/x in the second integral above, making a small rear- 
rangement and performing two elementary integrations, 
the final expression for the integral ( [f5| ) turns out to be 



Next, applying the identity 

1 \ T . fl 



Li 2 



As for the integral one can write it in the form 



G 2 (xi,X2) 



o a; — xi 



ln(x) 



Ax- 

X-(l-Xi) 



ln(. 



~*.\X) 

ax- 

o ^ - x 2 



ln(x) 



o dX x - (1 - x 2 ) 



, (50) 



which, when expressed in terms of the Euler dilogarithms, 
takes the form 



G 2 (xi,x 2 ) 



-Li 2 
Li 2 



+ Li 2 
- Li 2 



1 



1 — X\ 
1 

1 — X2 



(51) 



in (pil), we arrive at 



G 2 (x u x 2 ) = -2 Li 2 ( — ) + 2 Li 2 (— 

X\ I \x 2 



2 V -xi 



2 V -sea 



.(52) 



The remaining integral (E7h is elementary, and given by 



G 3 (xi,x 2 ) = In (xix 2 + iesign(/i|)) 



In 



1 — Xl 
-Xi 



In 



1 - x 2 
-x 2 



(53) 



Upon combining Eqs.(^J), Q4S|) , (|52j) and (^3|), we arrive 
at the final expression for the IR finite triangle scalar 
integral given by Eq. (|38|). 
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